SOME REMARKS ON SINC INTEGRALS AND THEIR CONNECTION 
WITH COMBINATORICS, GEOMETRY AND PROBABILITY 
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Abstract. We give an alternative, combinatorial/geometrical evaluation of a class of 
improper sine integrals studied by the Borweins. A probabilistic interpretation is also 
noted and used to shed light on a related combinatorial identity. 



1. Introduction 
In [1] , among other things the Borweins study integrals of the form 




where sinc(x) := (sinx)/x for x 7^ and sinc(O) := 1. They use a version of the Par- 
seval/Plancherel formula to prove that, subject to certain conditions on the parameters 
dj, the sequence I\, I2, ■ ■ ■ is decreasing; they also give an explicit evaluation of /„ by 
expanding the product of sines into a sum of cosines and then integrating by parts. As 
the Borweins observed, the integral (11 .ip can be interpreted geometrically as the volume 
of an n-dimensional polyhedron obtained by cutting an n-dimensional hypercube by two 
(n — l)-dimensional hyperplanes, and it is interesting to note that their formula is essen- 
tially a signed sum over the vertices of the hypercube. It may therefore be of interest to 
give an independent evaluation of (11.11) using methods from combinatorial geometry. The 
derivation is provided in £j3l 

In [2], related integrals arise in a probabilistic setting — namely the problem of de- 
termining the probability density of a sum of independent random variables uniformly 
distributed in different ranges. We employ this interpretation here to give an alternative 
proof of an elegant combinatorial identity related to the evaluation of I n using methods 
from probability theory. See Proposition [1] below. 
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2. Fourier Background and Connection with Probability Theory 

For our purposes, it is more convenient to take products from 1 to n, as opposed to 
to n in ( 11. II) . Thus, we fix a positive integer n, and let a = (ai, 02, ... , a„) be a vector of 
positive real numbers. For each j — 1, 2, . . . , n, define a step function %j : R — >• R by 

Clj, 

a v (2.1) 

a r 

Let /„ : R — > R denote the n-fold convolution ^i * X2 * • ■ ■ * Xn; so that for all real x, 

/oo /*oo /*oo 

Xi(^-ya) / X2 (2/2-2/3) ■ ■ Xn-i(y»-i-yn)xn{y n ) dy 2 ... dy n 
00 j —00 j —00 

It is now easy to see that f n (x) is simply the volume (in the sense of Lebesgue measure 
on R n_1 ) of the (n — l)-dimensional polyhedron 

{(xi, X2, ■ ■ ■ , x n ) e R n : YTj=i x j = x anc l \ x j\ < a j f° r j — 1) 2, . . . ,n} (2.2) 

divided by the volume 11^=1 °f the n-dimensional hypercube 

{(xi, x 2 , . . . , x n ) G R n : \xj\ < aj for j = 1,2,..., n}. (2.3) 

The volume of the region (12.21) will be evaluated directly in §21 but first we connect this 
problem with that of evaluating integrals such as (II. ip . The Fourier transform of f n is 
given by 

Ut) = Y[xj(t) = II / eltX Xj( x )dx = J]sinc(a^), t G R. (2.4) 

i=i i=i - 7 - 00 i=i 

Since for all real x, /„(x) = |/„(x+) + |/„(x— ) by continuity of f n for n > 1, and by 
definition of xi when n = 1, Fourier inversion gives 

1 /"OO 1 /»oo ^ 

fn(x) = — / e- tte / n (t) dt = — e Ux TT sinc(a^) rft, x G R, (2.5) 
2tt > /_ 00 2tt 

where we have made the change of variable t *—>■ —t and used the fact that sine is even. 
We remark that apart from our convention of starting products at j — 1, / n (0) coincides 
with the Borweins' integral (II. ip up to a numerical factor. More generally, f n (x) falls 
under the scope of their Theorem 2(ii). However, to facilitate the presentation of our 
combinatorial evaluation, it is convenient to reformulate their result as follows. We first 
make the following 

Definition 1. Let r : R — > R be given by 

1, if x > 0, 

r(x) = { |, if x = 0, (2.6) 
0, ifx<0, 
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and for x real and n a positive integer, let := x n ~ l r{x). Note that x Q + = t(x) and 
x™ = (max(i, 0))™ for n > 0. 

We are now ready to state our reformulation. 

Theorem 1. Let n be a positive integer, let a = (a 1; a2, . . . , a n ) be a vector of positive 
real numbers, and let f n {x) be as in 112.5]) . Then for all real x, 



fn(x) 



L ££{-1,1}' 



3=1 



[n 



- 1)! \{{2a 3 ) 

3=1 



(2.7) 



in which the sum is over all 2" vectors of signs e = (ei, 62, ■ ■ ■ , s n ) with Ej = ±1 for each 
j — 1, 2, . . . , n, and of course e ■ a denotes the inner product YTj=i £ j a j- 

To see more directly the connection with Theorem 2(h) of [T], simply replace the r 
function in Theorem [T] by the sign function using the relationship sign(x) = 2t(x) — 1 for 
x real. It then suffices to prove the identity 



E 



[x + e ■ a 



r 1 n 
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igR, 
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££{-1,1}" 



which is a special case of [fl Theorem 2(i)], there proved most elegantly using the method 
of generating functions. We conclude this section with an alternative proof of (I2.8P based 
on probabilistic considerations. For the reader's convenience, the result is restated as 

Proposition 1 (Theorem 2(i) of [T]). Let n be a positive integer and a = (ai, a 2 , . . . , a n ) 
a vector of positive real numbers. Then for all real x, we have 



0. 



J2 {x + e-aY Y[ej 

££{-1,1}" j=l 



n\2 n Y[< 
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if r = 0, 1,2,. 
if r = n. 



Proof. As in [2], we note that f n (x) represents the probability density of the sum of n 
independent random variables Xi, X 2 , . . . , X n with Xj uniformly distributed in the inter- 
val [—Oj, a,j] for each j — 1, 2, . . . , n. In light of Theorem [U the cumulative distribution 
F n (x) := Pr(]TJ = i Xj < x) is given by 



F n {x) 



fn(y) dy 



E 1 

££{-1,1}™ 



x + e ■ a) 



3=1 



Re, / n! 11(2 



(2.9) 



Let A n := Y^j=i a j- Since each Xj is uniformly distributed in [— a,j, aj], the sum YTj=i Xj 
must fall within the interval [—An, A n ], whence F n (x) = F n (oo) = Pr(^?=i ~Xj < A n ) = 1 
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for all x > A n . But if x > A n , then x + e ■ a > for each e G {—1, l} n , and so for such x 
we can drop the subscripted "+" from (12. 9p . It follows that 

n n n 

(x + e-a) n Y[ej = n\Y[{2 aj ) =n\2 n Y[ aj (2.10) 

£€{-1,1}" 3=1 j=± 3=1 

holds for all x > A n . Since the left hand side of (12.101) is a polynomial of degree at most 
n in x which is constant for all sufficiently large values of x, it must in fact be constant 
for all real x by the identity theorem. This proves the case r = n of the proposition. The 
other cases follow by repeatedly differentiating with respect to x. □ 

Remark 1. A related probabilistic interpretation of the integral ( ll.ip . which is essentially 
our f n (0), can be found in [H p. 81, Remarks l(b)(ii)]. 



3. Combinatorial Proof of Theorem Q] 



We will prove the formula (12. 7p of Theorem [T] by evaluating the volume of the polyhe- 
dron (I2.2p directly. Let us denote this volume by V n (x). Since the case n = 1 is trivial, 
we will assume n > 1. Let the vector a be as in the previous section, and let N be the 
set of vectors in R ri whose respective components exceed the components of —a. Thus, 
iV := {(xi, x 2 , ■ ■ ■ , x n ) G R n : Xj > —dj for j = 1, 2, ... , n}. To each subset A of N, 
we associate a non-negative function w(A) : R — > R defined by 

w(A)(x) := Vol{(xi,a; 2 , . . . , x n ) G A : YTj=x x o = 

where, as in ( 12.21) . the volume is taken in the sense of Lebesgue measure on R n_1 . Clearly 
w is a non-negative additive weight-function on the subsets of iV; for if A and B are 
disjoint subsets of N, then w(A U B) = w(A) + w(B), and w(A) > for all A C N. 
Therefore, if we let S := {1,2, ... ,n} and := {(xi, x%, . . . , x n ) G iV : x fc > a k } for 
k G S, then 



P| (N \ Rj) = {(x x , x 2 , . . ■ , x n ) G R n : -aj < xj < dj for j 



1,2, 



and hence in view of (I2.2p and the definition of w, we have for every x G R, 



V n (x) 



w 



f|(iV\^ 



ix 



n}, 



(3-1) 



since the underlying sets differ by a set of measure zero in R n 1 . The Inclusion- Exclusion 
Principle states that 



w 



■jeS 1 TCS j£T 

where the sum is over all subsets T of S, \T\ denotes the cardinality of T, and under the 
usual convention for intersections over empty sets, the term corresponding to T = {} is 
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w(N). Fortunately, the weights of the sets Rk and their intersections are relatively easy 
to compute. We make a linear change of variable 

n 

yj'-=^2 x k, J = 1,2, ... ,71. 

k=j 

Then yj = jjj + i+Xj for j = 1, 2, . . . , n — 1, y n = x n , and the Jacobian of the transformation 

is 1. Given T C 5, let b = b(T) = (pi, b%, ■ ■ ■ , &n) be the vector whose jth component bj is 
equal to a,j if j G T, and 6j = — cij otherwise. Then 

R(T) := f^Rj = {(x 1 ,x 2 ,...,x n ) G R n : x, > bj for j = 1, 2, . . . , n}. 

For each j — 1, 2, . . . , n, let Cj := X]fc=j V If £ £ ^C0> then ^ > Cj, and y^- — 6j > > 
Cj + i for j = 1,2, ... ,7i — 1. Thus, 

rVi—h pyz—b% rVn-i—bn-i 

w ( R ( T ))(yi) = / •••/ dy n - ■ -dy 3 dy 2 . 

Now make the change of variable Zj = yj — Cj for j = 1,2, ... ,n. Then each Zj > 0, and 
we obtain 

z-n-i „n— 1 



/■zi rz n -i 

w(R{T))(y 1 ) = / / ••■ / dz n ---dz 3 dz 2 
Jo Jo Jo 



(n-l)f 

In light of the constraint y\ — c\ = z\ > 0, we see that 

(yi-cor 1 _ (yi-E"=ifcj 

(n-l)l (n-1)! 
Thus, from (13 .ip and (13.21) it follows that 

1 n 

K(x)=^(-i)' T ^(n^)(^)=7^w e (^+^a)r i n^ 

where we have made the correspondence Ej = — 1 if j G T and Ej = 1 if j ^ T in each 
term of the sums. Since = V n (x) / Y\^ =l (2dj), the proof of Theorem CD is complete. 

□ 
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